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We study the time dependent potential energy W(t) — U(x(0)) — U(x(t)) of a particle diffusing 
in a one dimensional random force field (the Sinai model). Using the real space renormalization 
group method (RSRG), we obtain the exact large time limit of the probability distribution of the 
scaling variable w = W(t)/(Tlnt). This distribution exhibits a nonanalytic behaviour at w = 1. 
These results are extended to a small non-zero applied field. Using the constrained path integral 
method, we moreover compute the joint distribution of energy W(t) and position x(t) at time t. In 
presence of a reflecting boundary at the starting point, with possibly some drift in the + direction, 
the RSRG very simply yields the one time and aging two-time behavior of this joint probability. It 
exhibits differences in behaviour compared to the unbounded motion, such as analyticity. Relations 
with some magnetization distributions in the ID spin glass are discussed. 
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I. INTRODUCTION 



The description of glassy activated dynamics is one of the challenges of the present theory of glasses and disordered 
systems |]||. In fully connected, mean field approximations of these systems, the energy initially decreases and 
saturates at a threshold: the asymptotic large time dynamics then remains entirely dominated by high saddles in 
energy landscape and not by barriers. A recently studied model H, still mean field in character but which incorporates 
some effects of finite connectivity exhibits an additional regime of late times with subthreshold dynamics. It mimics 
an activated dynamics between lower and lower energy states, expected in any realistic physical system, leading to an 
interesting, slow but unbounded, decrease of the energy. 

At the other end of the spectrum, the Sinai model Q| of a particle diffusing in a one dimensional quenched random 
force field is dominated by activated dynamics and, despite its simplified character, is a non trivial example of a 
glass. It has the advantage of being analytically tractable, most notably through the real space renormalization group 
method (RSRG) and by various other methods Its extensions to many particles describe a variety 

of coarsening, domain wall growth and reaction diffusion models with disorder in one dimension, also analytically 
tractable via the RSRG [ ^6p7|JI^ , ^5| | . It is thus interesting to study, within the Sinai model, how the energy decreases 
as a function of time, as the particle jumps from deeper to deeper wells. 

The aim of this paper is to obtain exact results on the large time behaviour of the energy E(t) as a function of time of 
a particle in Sinai model. The energy lansdcape, denoted U (x), is itself a random walk with (U (x) — U(x')) 2 ~ \x — x'\ 
at large scales. Fixing the energy to be zero at initial time, and denoting x(t) the trajectory of a single particle, it is 
defined as: 



E(t) = U(x{t)) - U(x(0)) = -W(t) (I) 

i.e it is minus the loss of potential energy W(t), the quantity on which we focus from now on. In the Sinai model 
this quantity exhibits non trivial sample to sample fluctuations, even in the large time limit. In this paper, we will 
thus be interested in the distribution D t (W) of W over both thermal histories and environments (assuming a uniform 
probability to start at any site). To our knowledge, this probability distribution of the energy has not been computed 
previously. 

We will first consider the unbiased diffusion and use the RSRG approach (Section [[I B| ). This method consists in 
coarse-graining the energy landscape in a way that exactly preserves the long time dynamics lO.Rj. One decimates 



iteratively the smallest energy barrier in the system stopping when the Arrhcnius time to surmount the smallest 
remaining barrier T is of order the time scale of interest t, with V = Tint. The large time dynamics becomes 
asymptotically identical to the so called "effective dynamics" where the particle is with probability one at the bottom 
of the valley in the renormalized landscape containing the starting point. As we find the energy at this point is such 
that W(t) ~ Tint, i e the decrease is only logarithmic in time (much slower than the decay of the energy in the model 
studied in Ref ||). The correct scaling variable, which exhibits a universal probability distribution at large time is 
then: 

W{t) 

W= TnVt (2) 

and below, we compute this universal large time distribution D(w). The most notable result, besides the simple, 
exponential form of this distribution, is that this function is found to be non analytic at w = 1. For any finite time it 
is of course analytic, but in the large time limit the distribution of the rescaled quantity exhibits a jump in its second 
derivative. As discussed below this can be traced to the decimation process in the RSRG, but since this is a property 
which can be measured with no a priori knowledge of the physics of the system, e.g. via a numerical simulation, it 
provides a rather remarkable " experimental" signature that the glass fixed point describing this system is of the type 
exactly solvable by RSRG, i.e. characterized by an infinitely broad distribution of barriers. A interesting question is 
whether this is a more general property: this suggests to search for obscrvablcs exhibiting similar signatures in other 
systems with "infinite disorder" glass fixed point [l5| ]. 

The decrease with time of the energy can be related to the if-theorem for this system. One recalls that the 
generalized free energy: 



G{t)= I dx[TP(x, t) In P(x, t) + U(x)P(x, t)] (3) 
is always decreasing in each sample : 
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a simple consequence of the Fokker Planck equation for the probability distribution P(x, t) of the position of the 
particle dtP — —d x J, with J = —Td x P — Pd x ll. Here one can easily see, since the effective dynamics becomes 
asymptotically exact, that: 

G{t) „ 
lni t—>oo lnt 

confirming the considerations given in jl8j . 

Next, we will study the case of a small additional applied external force, for which the RSRG is still exact, and 



the result for D{W) will be generalized to that case in Section II C. In the following Section III we use a different 
and complementary method, the constrained path integrals, to obtain the joint distribution of the position x(t) and 
energy W(t). This yields information on the distribution of positions corresponding to a fixed potential energy loss, 
and we show, in particular that this distribution becomes deterministic at large W. 



Another situation studied here in Section IV corresponds to the geometry of the half-line, with a reflecting wall at 
the starting point. Then the RSRG allows to solve not only the joint distribution of position and energy at time t 
but also the full aging two-time behavior of this joint probability. The resulting physics is slightly different, however, 
as the distribution of energy does not exhibit the nonanalytic behaviour. 

Finally, the quantities studied here also have an interpretation in terms of a spin model, the ID spin glass. We 
show in Section [v| that the distributions of magnetizations of that model on some intervals is related to the above 
results. 

II. DISTRIBUTION OF ENERGY 

In this section, we compute the distribution D t (W) of the energy W(t) (Q) by solving the appropriate RSRG 
equation. 

A. Review of the main results from the real space renormalization method 

We now briefly recall the RSRG method for Sinai landscape, and we refer the reader to (ll]] for details. In the 
RSRG, the original energy landscape is replaced by a renormalized landscape at scale T, via an iterative decimation 
procedure where barriers smaller than T are eliminated. This generates a joint distribution Pr(F, I) of the barrier F 
and length I for the renormalized bonds in the landscape at scale T. 

In the symmetric case, the probability distribution of the rescaled variables rj = and A = ^ flows towards the 
fixed-point [f[l[ 

P*(fi. A) = LT7\ ( ^ e -v^coth^\ (6) 
v ; Vsinh Vs / 

In particular, the distribution of the barrier alone is a simple exponential 

r+oo 

P*( V )= / dAP*(7 ? ,A) = e-" (7) 
Jo 

In the effective dynamics |Tl|], the particle starting at x — at t — is at time t in a finite neighborhood |l4| ] 
around the bottom of the bond of the renormalized lanscape at scale 

r = Tlnt (8) 

which contains the starting point x — 0. After rescaling X = j^, the effective dynamics gives the asymptotic exact 
result of the Kesten distribution |ll| for the rescaled position (diffusion front). Similarly here, the effective dynamics 
will give the asymptotic exact result in the limit t — > oo for the rescaled variable of energy 

w w 
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B. Distribution of energy for the symmetric Sinai diffusion 



We introduce the probability Dr (F, W) that the point x = belongs at scale V to a bond of barrier F and is at 
potential W above the lower potential of the bond, i.e. U(x = 0) — U m i n = W. It is normalized as: 



rF -i roc 

/ dWD r (F,W) = = dllPr(F, 
Jo h Jo 



I) (10) 



since the probability that the origin x = belongs to a renormalized bond of barrier F and length I is lPp{F, l)/lr- 
Within the effective dynamics, we may compute the probability distribution D t {W) of the energy W at time t as 

D t (W) = / dFD v {F,W) (11) 

The consideration of the various effects of the decimations processes (see figure [j]) leads to the following RG equation 
for Dr(F, W) 

d r D r {F,W) = -2D r (F,W) [ Pr{TJ 2 ) (12) 

Ji 2 >o 

+ [ f P r (F 1 )P r (T)D r (F 3 ,W)S[F-(F 1 +F 3 -T)} 

+ [ [ [ DriF^W^PrirjPriF^diF-^+Fz-rfiSlW-iW' + F 3 -T)] 

JFi>T Jf 3 >F J0<W'<Fi 

+ [ f f PriF^DriT^^PriF^SiF-^+Fz-TMW-iW + F 3 -T)] 

Jf, >r J F,>r Jo<H"<r 




FIG. 1. Decimation of a bond of barrier F that gives a renormalized bond of barrier F — Fi + F3 — V : the point O can be 
either on the bond F3 that doesn't change the energy W = W', or on the bonds Fi or T that both lead to the incrementation 
W = W + F 3 - T of the energy. 

In the rescaled variables r\ = and w = -yr, the probability distribution T>r(r],w) = T 2 Dr(F, W) flows at large 
r towards the fixed point T>*(r), w) satisfying the equation 

0=(l + T))d v V*(ri,w)+wd w V*(ri,w)+ T drt'e'^-^V* (r)' ,w) (13) 

J rnax(0 7 w—l) 

+ dr ] 'e- T >'v*(ri-ri',w-ri') + e- r > dr]'V*(0, w - rf) (14) 

JO Jmax(0,w—1) 

in the domain 77 S [0, +00) and < w < 1+rj. Note that the probability that the origin belongs at scale T to a bond 
of barrier 77 is already known (|^) 

r 1+v i r°° 1-1-271 

/ dwT>*(r),w) = = d\\P*(r},\) = —— -e"" (15) 
Jo A Jo 3 
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It is now more convenient to introduce the variable v = 1 + r) — w instead to 77, and to set V*{rj,w) = e n Q{v — 
1 + r\ — w, w). The new function w) is now symmetric in (v, w) and satisfies 

= vd v $(v, w) + wd w $(v, w) - (v + w)§(v, w) (16) 
+ / dv'<S>(v',w) + dw'$(v,w') + dv'$(v',l-i/) (17) 

J max(0,l—w) J max(0,l — v) J max(0,l — w) 

on the domain {v > 0, w > 0, v + w > 1}. The function J(v, w) = d v d w <$>(v, w) satisfies the simpler equation 

Q=(vd v + wd w + 2-(v + w))J{v,w) (18) 

and thus takes the form 

v e v+w 

J(v,w) = C(-) (19) 

w vw 

where C(x) is an arbitrary function satisfying C(x) = C(-) to respect the symmetry (v,w). Since we are looking for 
a solution that is not exponentially growing at (v,w) — > 00, it is necessay to have C(x) = 0. The function <5>(w,w) is 
thus separable, and taking into account the symmetry in (v, w), we have 

$(v,w) = <p(v) + <p(w) (20) 

(21) 

where the function 4>(v) has to satisfy the equation obtained by replacing the form ( |2p| ) into ( |l6| ) 

= v<//(v)-(v + max(0,l-v))</>(v)+ dv'c/)(v')+ dv'<f>(v') (22) 

J maa:(0,l — v) Jmax(0,l — v) 

More explicitly, the equation becomes 

= V(j>'{v) ~ (f>(v) + 2 [ dv'(f>(v') for < v < 1 (23) 

= v<j>'{v) - v<j>{v) + [ dv'(f){v')+ ( dv'(j){v') for v > 1 (24) 
Jo Jo 

It is more convenient to derive these equations with respect to v 

= v<t>" (v) + 24>{v) + 20(1 - v) for < v < 1 (25) 
= v<j>" (v) + (1 - v)4> (v) for v > 1 (26) 

and to keep the boundary condition at v = 1 

O = 0'(l)-0(l) + 2 / dv'(f>(v') (27) 
Jo 

In the domain v > 1, the only solution that is not exponentially growing at v — ► 00 is the constant solution 

4>{v) = 0(1) for w > 1 (28) 

and the condition at v = 1 becomes 

0(1) = 2 / du'0(v') (29) 



The constant 0(1) has to be determined through the normalization condition ([L5j 



l + 2r? 
3 



I+77 rl+V 

dw$(v = 1 + r\ - w, w) = 2 / dw4>{w)=2l dw<p(w) + 20(l)ry (30) 



which leads to 
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(1) = I (31) 
1 1 

dv4>{v) = - (32) 

We know consider the equation in the domain v G [0, 1]. Since the equation is not local for <fi(v), it is convenient to 
introduce the sum 

S(v) = <f>(v) + <j)(l - v) (33) 

that satisfies the simpler equation 

v{l-v)S"{v) + 2S(v) = (34) 
Two linearly independant solutions of this equation are 

S sym (v) =v(l-v) (35) 

g-antfeym/^ = 2v (\ - v ) m + 2« - 1 (36) 

1 — V 

Since S{v) has to be symmetric in (i>, 1 — v) by definition, and has to be normalized to 

[ dvS{v) = 2 [ dvcf>(v) = i (37) 
Jo Jo 3 

we get 



and this corresponds to 



We thus finally get 



S{v) = 2v(l - v) (38) 



^) = i-^(l-«) 3 (39) 



4>{v) = - - -(1 - v) 3 for0<u<l (40) 
3 3 

<j>(v) = - for v > 1 (41) 
3 



The final result of this section is that, at large T, the probability 2?* (77, w) that the point x = belongs to a bond of 
^1 
r 



rescaled barrier r/ = F „ T such that the point x = has a rescaled potential w = J 7 """ above the lower potential 



of the bond, reads 

D, (77, w ) = e _,) $(l + 77 - w, w) = e~ v (<t>(w) + <f>(l + r\ - w)) (42) 
In particular, the distribution of the rescaled energy w IM) alone reads 



X>*H = / dr]V^ lW ) (43) 

J max(Q,w — l) 

= 6{w < 1) (4 - 2w - Ae~ w ) + 6{w > 1) (2e - 4) e~ w (44) 

This function is represented on Figure (^|). It is continuous at w = 1, as well as its first derivative, but the second 
derivative is discontinuous at w = 1 . The contribution to the normalization of the domain w > 1 reads 

r+co a 

/ dwVJw) = 2-- = 0.52848 (45) 
Ji e 

The mean value is simply 



Wit] [ + °° 4 
lim —4 = / dwwV*(w) = - (46) 
t^oo 1 ait Jn 3 
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FIG. 2. Exact asymptotic probability distribution D,(w) of the rescaled energy 



C. Distribution of energy for the biased case 



As discussed in details in Ref. fill , the RSRG method can be applied to the Sinai diffusion in the presence of a 
small external bias 5, and the results are valid in the double limit S — -> and T — > oo with the scaling variable 



7 = ST = STlnt 



(47) 



being fixed. 

The above study of the symmetric case is generalized to the biased case in Appendix [A|. The final result for the 
probability distribution D t (W) of W(t) is the sum of the contributions of the two types of bonds: 



D t (W) = D+(W) + D~(W) 



with T — Tint and 

Dp(W > T) 
for W > T, and 



d±(w < r) 



sinli <Sr((Scoth(5rT3(5) 



^ _ 2 e -r(=F<5+<5coth5r)\ e -(lV-r)(=F<5+<5coth<5r) 



2sinli Jr(coth <5r=F3) 



(-4e- W+a- cothsr) + 3 ^ CQth6T + (1 ± coth6T)e T2SW 



(48) 



(49) 



(50) 



for < W < r. 



As in the symmetric case, the probability distribution D^iW) and its first derivative dwD^t (W) are continuous at 



W = T, but its second derivative d^D^ (W) is discontinuous at W = T. 
The full normalizations 



dWD±(W) 



5± 2 7 _ 1 zp 2 7 

4 sinh 2 7 



(51) 



correspond as it should to the ratios lp / (/jt + l r ) in terms of the mean lengths of the renormalized bonds. 
The mean value for the + case reads 
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r+oc 

W+(t)= / dWWD±(W) (52) 
Jo 

e 7 (l- 2e 7 - 7COth7 )(e 7 + ( (2e6T + ^^g^ ~ 1+27+472 } - 16 (7 + e 7 sinh 7 )e^ 3 -™ th ^) 



%oth7-3) &6(e 2 i-2) 



and varies between W + {t) — > (2/3)T as <5 — > and 



W + (t) ~ -e^ 7 (54) 
7^00 2o 



III. JOINT DISTRIBUTION OF ENERGY AND POSITION 

In this section, we compute the joint distribution of position x = x(t) — x(Q) and energy W — U(x(0)) — U(x(t)) 
by the direct path-integral approach introduced in Jl5| . 

A. Recall of the path-integral approach 



As explained in details in |15|, the properties of the renormalized landscape can either be obtained as solutions of 



RSRG equations or by the computation of constrained path-integrals. In this second approach, the basic path integral 



rVi=V rx -, ,y 

F [VatVb] (V ,x ,V,x)= DV(y)exp(- dy(-(—) 2 )er Va ] ({V}) (55) 

Jv m() =v Jo 4 <*y 

over the Brownian paths going from V(xq) — Vq to V(x) = V without having touched the boundaries V — V a and 
V = Vb can be expressed in Laplace transform with respect to [x — x$) as 



F[v a ,v b] (v, P \v ) = w ,.y a Vb \ *-<y>n v*)*+(Vo,r, H) if v a <v<v 



(56) 



F[v a ,v b )(V,p\V ) = 1 $ + (V,p;V b )$-(V ,p;V a ) if V < V < V b (57) 
W{p;V a ,V b ) 

where $_ and $ + are solutions of the differential equation 

dl$= P F (58) 

that vanish respectivly at V — V a and at V = Vb 

$-(V,p;V a )=sinh^p(V-V a ) (59) 
$+(V,p; V b )= wahy/p(V b - V) (60) 

(61) 

and where W{p\ V a , Vb) is their wronskian 

W(p;V a ,V b )= y/psmhy/p(V b -V a ) (62) 
As explained in jL5|, the probability distribution Pr(F, I) corresponds to the path integral 

B r (F,l) = jf^ DVex P (- J dx{-(—f)Q^{{V}) (63) 

where F > T and where 0p ({V}) constraints the paths to remain in the interval ]0, F[ for < x < I and to not perform 
any returns of more than T, i.e. for two arbitrary < X\ < x 2 < I, the potential has to satisfy V(xi) — V(x 2 ) < T. 
As shown in pi , it can be obtained from the path-integral ( |57| ) as 
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B T (F,p)= f* e -(F-r)VPcothVpr (64) 
sinh ^/pT 

in agreement with the fixed-point solution (Q) of the RSRG equations. 

Similarly, the probability Er(V,l,Vo) to go from (Vb,0) to (V,l) without making any return of more than T and 
where V > Vq is the maximum reads in Laplace [ fl5[ 

^r(X / ,p,y )=e- (y - Vo) v^ coth v^ r (65) 
We are now in position to compute the joint distribution of energy and position via a path-integral decomposition. 



B. Joint distribution of energy and position for the symmetric case 

The law is symmetric in x — > —x so we will consider x > and write the path-integral representation 

P t (x > 0, W)= ±= < 6(V(x) ~ W) > {vm (66) 

-, r+oo r+oo ,V(l) = F { 

^ / dF dl 2?V(y)e~* f ° dy ^> 5{V(x) -W)Q^{V} (67) 

Jx Jv(0)=0 



r 2 



where again Qr{V(y)} means that the path is constrained to remain in (0, F) with no descending segment of more 
than r. Since this constraint doesn't factorize into the two-sides y < x and y > x, we have to keep track of the 
maximum V m reached by V(y) for < y < x to impose that the path doesn't go below V m — T for x < y < I. 
We now compute the Laplace transform 

r+oo 

P t (p,W)= / dxe-**P t (x,W) (68) 



(69) 

for the two cases W > T and < W < F. 



1. Case W > T 

Let us introduce the following notations (see Figure ||) : We note V m the maximum reached by V(y) for < y < x, 
and x m the point where it is reached V(x m ) = V m . We note x' m the first time in x < y < I where V(y) reaches 
again V m . For the given parameters (y m ,x m ,x' m ) 1 we have to sum over the Brownian paths satisfying the following 
conditions : (i) they go from (0, 0) to {V m , x m ) without making any return of more than T, and V rn is reached for the 
first time at x m \ (ii) they go from (V rn ,x rn ) to (W,x) without touching the absorbing boundaries at V m — T and at 
V m + e. (iii) they go from (W, x) to (V m ,x' m ) without touching the absorbing boundaries at V m — T and at V rn + e. 
(iv) they go from (V m , x' m ) to (F, I) without making any return of more than T and F is the maximum. 

The summation over the Brownian paths satisfying all these constraints leads to 

^ pW+T fx 

P t (x,W>r)=-f dV m dx m B r (V m ,x m )(d Vo F lVm ^ r y m] (W,x\V ,x m ))\ Vo=Vm (70) 

1 Jw Jo 

dx' m {d Vl F [Vm ^ r y m] {V u x' m \W,x)\v^v m dF dlE r (F,l;V m> x' m ) (71) 

•>V m Jx' 



where the functions Fty at v b ](V,x\Vo,Xo),Br(F,l) and Er(F, I; V, x) have been defined in section (HI A) and are given 
by the explicit expressions ( |57| ) ( |64| ) ( |65| ) in Laplace transform. 
As a consequence, in Laplace transform with respect to x, we get 

P t (p,w> r)= ' f w+ \ Vm _^ e -^-r^^ ^ vp ( r (v - w)) 



r 



v„, 



sinhygr {w _ r)V p COth 

r 2 ^ 



(l - 2 cosh Vpr e - r v / P coth Vpr^j ( 74 ) 
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y 



FIG. 3. Properties of the renormalized bond for W > I\ 



Rescaling the variables 

r W 

we get the Laplace transform of the exact asymptotic distribution of the rescaled variables X and w 

r+oo 

V(s, w>l)= / dXe- sX V(X, w>l) 
Jo 

sinh \/s 



(eV^coth^ _2coshVi) g-^v^coth^ 



(75) 

(76) 
(77) 



In particular, the distribution of rescaled energy w alone reads (by taking into account the two sides X > and 
X < 0) 



V(w > 1)= 2V(s = 0, w > 1) = 2(e - 2)e 
in agrement with the RSRG result (E3). 



(78) 



Case W < r 



For VF < r, there are two contributions coming from the cases V m > T and V m < T, where V m is again the 
maximum reached by the potential V(y) for < y < x 



p t {x, w < r) = p t (1) (x, w < r) + p t (Z> {x, w < r) 



(2), 



(79) 



For V^j > T, we have the same decomposition as in ([7l|), except that here we have to integrate over V m € (r, W+T) 
(instead of V m e (W, W + T)) 



1 



w+r r x 



P t 1 >(x, W)= ^2 J dV m J dx m B T {V rn ,x m ) (dv F[v m -r,v m ]{W,x\V ,x m )) \v =v m 

r~\-oo p-\-oo p-\-oo 

dx' m {d Vl F [Vm _ T y m] {V llX ' m \W 1 x)\ Vl=Vm dF dlE r (F,l;V m ,x' m ) 



V m Jx 



(80) 
(81) 
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In Laplace transform, we thus get 



1 



p«(p,^<r)= — 



w+r 



dVr, 



Vp c -(v m -n,/pcoth,/pr sinh Vp( t - (}jn - w)) 

' sinh y/pT sinh ^JpT 



+00 



dFe 



(r - (v m - w)) 
r 

= J_ / smhy/ppT-W) _ w sinh ^p(T-W) _ sinh 2^pT & _ w ^ coth ^ 
r 2 V VP sinhy/pT y/p 



Rescaling the variables as in (|75[), the first contribution in (79) reads 

V {1) (s,w <l)= [ + dXe- sX V {1) {X,w < 1) 
Jo 



(82) 
(83) 
(84) 

(85) 



= sinh ^(2- w) _ ^ sinhyi(I - w) _ sinh2yi e _ Wv ^ cot hVi 
y^i sinh -y/s y^i 

We now consider the second contribution corresponding to the cases where the maximum V m reached by V(y) for 
< y < x satisfies V m < T (see Figure ||). We then have to sum over the Brownian paths satisfying the following 
conditions : (i) they go from (0, 0) to (W, x) without touching the bondaries at V = and V = T. (ii) they go from 
(W, x) to (r, z) without touching the absorbing boundaries at V — and at 7 + e. (iv) they go from (r, z) to (F, I) 
without making any return of more than T and F is the maximum. 




z I 

FIG. 4. Properties of the renormalized bond for V m < Y. 
The summation over the Brownian paths satisfying all these constraints leads to 



y 



P?\x,W)= (d Vo F] ,r] (Vb , 0; W, x)) Vq =q I dz(-dvrF [0in (W,x;,V',z)) vl= 



1 



F 2 

/+00 /<+oo 
dF J dlE r (T,z;F,l) 

and thus in Laplace with respect to x, we get 

1 f +oc 

p {2) (p,w)= -j (d Vo F [0 , r] (v ,w ]P )) Vo=0 (-d v ,F [QT] (w,v';0)) vl=r / dFE r (r,F ; o) 



(87) 



w 



sinh(r - W)y/p 
r 2 sinh Ty/p 



(89) 
(90) 
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Rescaling the variables as in (|75[), the second contribution in fl79p reads 

V^( S ,w<l)^ [ + °° dXe- sX V^(X,w < 1) = w sinh (l-^)^ (91) 
' smh Vs 







Adding the two contributions (86) and (Dlj), we finally get the exact asymptotic distribution of the joined rescaled 



variables (X, w) in the domain w < 1 

V{s, w < 1)= V (1) (s, w < 1) + 7? (2) (s, w < 1) (92) 
_ sinh-y/s(2- w) sinh2yg e _ wVgcotllA/i 

In particular, the distribution of rescaled energy alone w reads 

V(w < 1)= 2V(s = 0,w < 1) (94) 
= 2(2 - iu - 2e~ w ) (95) 

in agrement with the RSRG solution (fi"4|). 

3. Kesten distribution for the rescaled position alone 

A check of the above result is that one recovers the Kesten distribution, obtained in [jllj via the RSRG, upon 
integration over the energy. 

For the domain w > 1, the distribution of the rescaled position X has for Laplace transform (R6h 



r>+oo 

V>(s)= I dwP(s lW > 1) (96) 

, / v / i — — \ 

(97) 



_ sinh y/j / 1 _ 2e _vScothV^ 
s \ cosh y^S J 

whereas it is given in the domain w < 1 by ( j9l| ) 



7><(s)= / dwV(s,w < 1) (98) 
= - f 1 - cosh + 2 sinh 2 v^e-^ coth ^) (99) 



Adding these two contributions, we recover as it should the Kesten distribution |fl| 



^. Displacements for fixed energy 

It is interesting to use the above solution to compute the moments of the distance traveled by the particle given 
that the (rescaled) potential energy lost is a fixed w = W/T\nt. The probability distribution of the rescaled distance 
X for a given rescaled energy q w (X) = V{X,w)J JdXV(X,w) has for Laplace transform @||) 

qw{s)= n*,u>l) = smh^i ( e ^coth^_ 2cosh ^ eW (i-VicothVi) for ^ > 1 (i l) 

~. . sinh y/s(2-w) sinha^/i -wJlcothJl 

V(S,W<1) e 

w (s)= - = — 7^— for w < 1 (102) 

W P(0, W <1) (2-«;)-2e-™ V ; 

The moments 
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c„M = | o + °° dXX n q w (X) =< >iu (103) 

represent the moments of the distance traveled over the set of trajectories and environments such that the potential 
energy lost is w. 

In particular, the first moment reads 

. , w 8 — 3e . , r«rn\ 

C1 W=3 + 6M) i0TW>1 (104) 

v ; 6(e<"(2- w) -2) v y 

and order n similar polynomial extensions are obtained for the n-th moment. As a comparison, note that c% is 5/12 
for the full Kesten diffusion front. 

These moments c n (w) are increasing function of w as expected since roughly speaking larger w correspond to larger 
displacements. For w = 0, they simply vanish c n (w) = 0, which correspond to 

g w =o(«) = l i.e. qw=o(X) = S(X) (106) 

For large w, the asymptotic behavior of the moments is simply given by 

c „h ~ (^y (to?) 



This interesting property can be understood by considering the large w limit for the distribution (101), which is 
dominated by the region s — ► : 



— >OG 



g lu ( S ) ^_e- w 3 (108) 

corresponding after Laplace inversion to 



w s 



q w {X) ~ (109) 

it; — >oo o 



i.e as it) becomes large, in terms of the original variables, we have 

x{t) _ 1 W(t) 
(Tint) 2 ~ 3 Tint 

with probability one. 



(110) 



C. Joint distribution of energy and position for the biased case 

The above calculation can be generalized in the presence of a small bias. 

The joint distribution of energy and position is given for the biased case by the following path-integrals with x > 

P t (±x,W)= =J-= < S(V±(x) - W) > {v±(y)} (111) 
t+ + 1— 

dF dl W^e-^o^wT 25 ) e r {V(y)}5(V(x)-W) (112) 

W Jx JV(0)=0 



sinh^ ST 



The explicit computation is done in Appendix ( [A 3| ) and leads to the following final results in Laplace transform 
for W > T and W < T respectivly 



V^L e ^ sr - 2( T S + ScothST)e- r ^ s+ScothSr '> 

P?(p, W > r)= \ si " h V^r ^ = ^ = e -(w-r )(T « + ^T^coth V^r) (n3) 

sinh 2 <5r ( +4^ 2 T45v^T^cothv7+^r) 

v sinh 2 r v ' 
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P?(p,W<T) = 



sinh ST 



p+s 2 



lh 2 y/p+S 2 T 



4S(+8 + y/p + S 2 coth y/p + S 2 T) 



[-2(+S + V / P+^ 2 coth ^T^r)e- M/(T5+ VP+^ coth Vp+^ r ) 

+2(^8 + \Jp + S 2 coth y/p + S 2 T) cosh y/p + S 2 W cosh SW 

-{y/p + 6 2 {l + coth 2 VP+^r) T 2(5 coth VpT^r)) sinh Vp + ^WcosIkW 



+ t(p ^ (T<5 + + 3 2 coth Vp + ^r) 2 ) cosh y/p + 8 2 W sinh SW 
a 



-( T S + y/p + 5 2 coth y/p + S 2 T)(~2 coth y/p + S 2 T - 



y/p + S 2 



8 sinh 2 VP+^r 



) sinh y/p + PW sinh <W] 



<5 sinh SW sinh(r - W) y/p + 5 2 
sinh 2 £T sinhT^p + 5 2 



(114) 

(115) 
(116) 
(117) 

(118) 
(119) 

(120) 



For the special case p = 0, these expressions coincide as it should with the results (|49j) and ([50|) corresponding to 
the distribution of the energy W alone. 

On the other hand, the integration over the energy W gives 



dWP^{p,W> T) 



y/p+S 2 



sinh T 



e T<5r _ 2 ( T s + S coth ST)e- r ^ s+s coth 5r > 



sinh 2 ST (- 



a 2 - 



lh 2 y/p+8 2 T 



E==- + AS 2 =f 46 y/p + 6 2 coth y/p + S 2 T)(^6 + y/p + S 2 coth >/p + (5 2 r) 



and 



H <± (P)= 



+r 



(/^(p,^ < r) 



5 2 



psinh 2 5r(( x /p + <5 2 coth y/p + 5 2 T - 25) 2 -p- 8 2 )) 



[2pe 



-r(T$+y/p+$ 2 coth y/ p +S 2 r) 



y'p + <5 2 sinh (5 



(£ 2 + 3<5 2 coth ST - 8 y/p + 5 2 coth Vp+^ 2 r(3 + coth tfr) 



P + S 2 



S sinh a/p + 8 2 T 

S 2 f _ y/p~+~8 2 sinYi8T \ 
p sinh 2 ST \ 8 sinh y/p + 8 2 T J 

whose sum coincides as it should with the distribution of the position alone [fTi"| 

P t ± (p)=P> ± (p) + P< ± (p) 

S 2 



sinh 2 y/p + 8 2 T 



)] 



p sinh ST 



Vp+ s2 c tat 



1 



sinh y/p+8 2 T 



\ 



(±8 + y/p + S 2 COth y/p + PT) 



J 



(121) 
(122) 



(123) 
(124) 

(125) 

(126) 

(127) 
(128) 



IV. RESULTS IN THE PRESENCE OF A REFLECTING BOUNDARY 

In this section, we consider the case where there is a reflecting boundary at the starting point x — jnj, with 
possibly some bias. It is interesting as it can be studied easily using the RSRG. 
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A. Joint distribution of energy and position on the half-line 



In the case of a bias, the joint probability distribution of E^(W, I) of the energy-loss W = U(t = 0) — U (t) and the 



distance I = x(t) — x(0) satisfies the following RG equation 

/>oo 

drE±(W,l) = -E±(WJ) / dl'P^(T,l') 

Jo 

pW />oo poo />oo />oo 

+ / dWt / d/i / dZ 2 / dF 3 / dl 3 E±(W 1 ,h)PT(r,l 2 )P±(F 3 ,l 3 ) 
Jo Jo Jo Jr Jo 



5{W - (Wi + F 3 - r))*(J - (Zi + z 2 + z 3 )) 



(129) 

(130) 
(131) 



For large T using the properties of the fixed point solution (|Al[ ) for P ± , and the properties (A4) of the functions 
U and u, this can be rewritten in the Laplace variable with respect to the length as 



drE±(W,p) = -uT(0)E±(W,p) + U+{ P )Uf( P ) / dW 1 E^(Wi,p)e-^ w - Wl ^^ p) 

Jo 

and we finally get that the joint distribution of position I and energy W takes the very simple form 
The distribution of the energy W alone is a simple exponential 



(132) 



(133) 



(134) 



contrary to the case studied in previous sections where the diffusion takes place on the full line. The absence of 
nonanalytic behaviour is traced to the fact that here decimations of the bond near the boundary do not occur. 



B. Aging for the joint distribution on the half-line 

It is now possible to obtain the correlation of the energy at two different times, i.e. the aging of the system in 
energy. 

In the case of a reflecting boundary at the origin, the single time diffusion front is simply 



P(W, x', t'\0, 0, 0) = E+ (W, x) 



(135) 



whose Laplace transform is given in equation (133). Similarly, the two-time diffusion front P(W, x, t, W , x', t'\00) = 
Ep r ,(W,x,W',x'), satisfies the RG equation given in (131) but with the initial condition at T = V : 



E r=r> v ( W > x ' W '> x ') = S ( W ~ W)S(x - x')E+, (W\ x') 



(136) 



Within the effective dynamics, this system has the flavor of a directed model, since W — W and x(t) — x(t') are 
always positive. It is thus convenient to define the Laplace transforms 



E+(X,p, W',p' 



dx'e- p ' x ' I dxe- p{x - x ' } / dWe- X(w - w '^E+ r ,(W,x,W',x') 

Jw 



(137) 



Using the fixed point solution (Al) for Pp together with the properties (A4), the above RG equation simplifies into 



9 r lnP+ r ,(A,p, W',p') = -%(0) + ElMErM = n v 1„ 



X + u+ip) 



4(Q) 

X + u+{p) 



with the initial condition 



E+ =T r ,(X,p, W',p') = E+,(W', P ') = u+(0)e 



(138) 



(139) 



We thus obtain 
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E+ TI {\,p,W\p')=E+(W\p')^ ^^TT ,NlM 



or more explicitly in terms of (W, W) and (p,p') 

E£ ir ,(W,p,W',rf) = u+{0)e- w ' u ti (p,) (6(W-W) + (u+,(p) - u+(p)) e - (w ~ w "> u r(p)^ (141) 
In particular for p — = p', we obtain the aging properties of the energy alone 

E+ r ,(W,W) = u+(0)e- w ' u t'W (s(W-W) + (u+(0) - u+(0))e" (w/ - w ' ) "i t(o) ) (142) 
For the symmetric case (5 = 0), it simplifies into 

Err(W,W) = ie-^'/r' ( S (W-W>) + (I - i)^"^')/^ (143) 



(144) 



r v v T' r 

In particular, the two-time Energy Correlation reads 

r+°° r+°° 1/1 1 \ 

U(t)U(t>) = dW dW WW'E+ r ,(W,W) = — - — - + — - 
Jo Jw> uf,(0) \Up/(0) (0)7 



In the limit 7 = T5 = T5\nt > 1 7' = T'S = T5\nt' 3> 1, the influence of the wall vanishes in that limit, and the 
result should coincide with the result of the full model without a wall. 

V. DISTRIBUTION OF THE MAGNETIZATION IN THE ID SPIN GLASS 

Finally, let us indicate how the above results can be translated in the language of the one dimensional Ising spin 
glass in an external field. We refer the reader to [|l6| for the details of the study of the spin glass by the RSRG 
method. We consider the Hamiltonian 

i=JV-l i=N 

H = - JiWi+i - ha t (145) 

with a bimodal distribution of random bonds Ji = Jei . In the absence of a field h = the spin glass has two ground 
states ±<7° with of = ei-.e<_i. In the presence of a field h the ground state of the spin glass is made of domains of 
either zero field ground states. At finite time t after the quench, these domain grow until they reach their equilibrium 
size at equilibration time T eq = Tin t eq = 4 J. Each of these domain consists of an interval between two frustrated 
bonds (i.e. such that JjUjCTj+i < 0). During the coarsening process the number of such bonds decrease with time. 

Choosing a point O randomly in space one can thus, at time t, corresponding to the RSRG scale T = Tint denote 
by i = R (and i = L respectively) the closest frustrated bond to the right (respectively the left) of point 0. In the 
spin glass the (absolute value of the) magnetization plays the role of the energy barriers in the Sinai model, more 
precisely: 

|M i0 | = | J2 **l = l E °i\ = -5j:\V(L)-V(p)\ (146) 

L<i<0 L<i<0 

where V is a Sinai type potential. It is then easy to see that if one defines 

\M L0 (t)\ \M 0R (t)\ 

Wl = -f^r W2 = ~T^r (147) 

then the joint distribution P(wi, W2) of the scaling variables u>i, W2 has been determined in Section II: 

P(w u w 2 ) =V«{r, = Wl +w 2 - 1, Wl ) = e- u,1 - u,2+1 (0(u)i) + <j)(w 2 ))9( Wl +W2-I) (148) 

where <fi(w) is given in (flT|). It exhibits, in particular, the same nonanalyticity as the distribution of the potential 
energy loss of a particle in the Sinai model. The distribution of the total magnetization |Mi,o(i) + Mqr(£)\ = 
\Mw(t)\ + |-^ofi(i)| is just the barrier distribution, i.e. a simple exponential, as found in |lq]. Note also, that the 



joint distribution of the equilibrium magnetizations is also (148) with the definitions (147) replacing Tint — > 4 J 



16 



APPENDIX A: GENERALIZATION TO THE BIASED CASE 



1. Recall of main results of the real space renormalization method 

In the presence of an external drift, there are two probability distributions PfT (F, I) describing respectivly the 
ascending and the descending bonds. The asymptotic distributions read [Tlfl 

Pr(F,p) = C/ r t (p)e- (F - r) "r(p) (Al) 
(p) = Vp + S 2 coth \Ty/p + S 2 } t S (A2) 



It is also usefull to recall the evolution equations 

d r 4(p) = -U+(p)Ur(p) (A4) 
drU±(p) = -i4(p)U±(p) (A5) 

In particular, the distributions of the barriers alone are in terms of the variable z — F — T 

P±(z)=u±{0)e~ zu rW ( A6 ) 

Finally, it is also convenient to introduce the density of renormalized valleys 

n r = _L= = , f = u+ (0)«F (0) (A8) 
1+ + /p smh (Td) 



2. Distribution of energy via the RSRG equation 

We introduce the probability Dp (F, W) that the point x = belongs at scale T to a ± bond of barrier F and is at 
potential W above the lower potential of the bond (i.e. U(x = 0) — U m in = W). They are normalized to 



pF -i />oo 

/ dWDy (F, W) — = — = / dllP±(F, 

Jo li + i r Jo 



I) (A9) 



The RG equations for Df {F, W) are 

<9 r £>r ( F , w ) = - 2D r ( F , w ) p r ( r ) ( A1 °) 

pOO 



/oo poo 
dF 1 dF 3 P^{F 1 )P^{T)D^(F 3 ,W)6[F - (F 1 + F 3 - V)] (All) 

poo poo pF\ 

/ dF ± dF 3 dW'D^F^W^i^P^iF^SlF-iF.+Fs-rmW -(W + F 3 -T)} (A12) 
Jr Jr Jo 
/>oo poo />r 

/ dF 1 / dF 3 / dW'P^iF^E^iT, W')P±(F 3 )5[F - (Fj + F 3 - T)}S[W - (W + F 3 - T)} (A13) 
Jr Jr Jo 



Using the variable z — F — T and the fixed point solutions Pf (z) (A7), the RG equations become 



(ft- - d z )D±(z,W) = -2u*D±(z,W) +itpu+ / dz'e~ (2 - 2 ' )u rD r (z', W) (A14) 

^maa;(0,lV-r) 

/•mm(.s,W^) />mm(z, W) 

+ufu+ / dz'e- 2 ' u r J D±( z - z ' ; VF-z') + (u r : ) 2 e~ zu r / dz'-Dff (0, W - z') (A15) 
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Using ( |A8| ) it is convenient to set 

D±{z,W) =nru^e- zu r->S!±(V = z + T -W,W) (A16) 
and to write the RG equations for tyi(V, W) 

[dr + (V + W -T)n T ]^(V, W) (A17) 

(i-V pW pmin(T.V) \ 

/ dV'V r z (V,W) + dW'V r z (V,W) + dV'*^(V',T-V')) (A18) 

Jmax(0,T-W) Jmax(0,T-V) J max{0.T-W) I 

Applying dydw to these equations lead to the simple equations 

[Br + (V + W - r)n r ] d v d w ^ (V, W) = (A19) 

Since we are loking for solutions that are not exponentially growing as (V, W) — » +00 , we are led, as in the symmetric 
case, to 

dvd w ^(V,W) =0 (A20) 

It is now natural to set 

^+(V,W) = #f (V, W) = $>r(V,W) = tpr{V) + ip r {W) (A21) 
and the equation for ipr(V) reads 

[Br + (V - T + maiz;(0, T - V))n v ] ipr{V) (A22) 

(/■V /.mm(I\V) \ 

/ dV'ipr(V) + / dV'MV) (A23) 

./ma:r(0,r-V) Jmai(OT-V) / 

We now differentiate this equation with respect to V 

[d r + {V -T)n r ]d v ^r{V) = for V > V (A24) 
BrB v ipr{V) = 2n r (#(V) +#(r- y)) for < V < r (A25) 

and keep the following condition at V = T 

dr#0Ok=r = 2n r / dVtyr(V') (A26) 
Jo 

The normalisation condition ( |A9| ) now reads 

£ dWD±(F, W) = n r J™ dllP^F, I) = -n T d p (uf^e^-^^ | p=0 (A27) 
= nr^(0)e-^- r )"r(o) {-d p \nU±(p)\ p = + (F - T)d P 4 (p)| p=0 ) (A28) 
Taking into account the changes of variables 

,F ,F 

/ dWD^(F,W) =nr^(0)e- (F - r)u r(°) / dWf±{F - W, W) (A29) 
Jo Jo 

+ r F 

= 2n r ur(°) e_(F ~ r)Ur(0) / (A30) 

Jo 

we get for F > T 

dW#(W) - -9plnf/±(p)| p=0 + (F — T)B p ]nu±(p)\ p = (A31) 
Deriving with respect to F, we get for any F > T 
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1 ± (sinh(2gT) - 2Jr ) 

2 p rWlp 85sinh 2 (<5r) 



For the value f = Twe get 



"dW^WO = -iaplnt^Wl^o = ^^(f) 1 (A33) 



The equation QA25D for 1^ > T is thus satisfied, as well as the condition (A26) at V = T. We now have to find the 



ipr(V =T) = Vt(T) given in (A32) and d v ip v {V)\v=r = 0. It reads 



solution ipr(V) to the equation (A25) on the domain < V < T that satisfies the continuity boundary conditions 



id sinn (ol J 

The final result is thus that the joint distributions Dfr (F, W) are given by 

£>±(F, = n r wre~ (i? ~ r)u r + Vr^ - WO) (A35) 

where 

^r(W) = -/ir(r) - d(T - W)h r (T - W) (A36) 

in terms of the auxiliary function 

(sinh(2<5VO - 25V) 
4<5sinh 2 (5r) 

We now compute the distributions of W alone 

r+co 

D±(W)= dFD±(F,W) (A38) 

Jmax(r,iv) 

/■min(r,wo 

= n r (/i r (r) - 0(r - W^)/i r (r - W 7 )) e" max(o,w-r> r _ np / dVujW 1 '-^'* ft r (y) (A39) 



My) = ^y:]^::** (as?) 







More explicitly, for W > T, we have a simple exponential form in W 



whereas for < W < I\ we have 



ftr(r)-u^ y dFe- (r - y)u ^r(^) j (A40) 



D r z (W<T) = n r ^h r (T)-hr(T-W)-u^ dVe (v ~ w)u r h T {V)^ (A41) 



These expressions lead to the final results (f49h and (pO) given in the text 



3. Joint distribution of position and energy 

The path-integral approach to compute the joint distribution of position and energy can be generalized to the biased 
case as follows. 

The basic path-integral becomes 



V,=V rx 1 d y 



\± atVb] (V ,x ,V,x)= _ DV(y)exp(- dy{-{— T 28Y)6^ Vb] {{V}) (A42) 
= e- 52 (*- x °W y - v ^F [Va>Vb] (V , x , V, x) (A43) 
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and thus in Laplace transform with respect to (x — xq), we have in terms of (|57 

F± a<Vi] (V,p\V ) =e ±s ( v - v °)F [Va>Vb] {V,p + 5 2 \V ) (A44) 
Similarly, and (p5h become 



B±(F v)= c ±sr VP + ^ -(F-r)( T f+V^coth^4^r) (M5) 

r sinh^pT^r 

^r(^p|F )= e -(v-Vo)( T 5+V^coth VS^r) (A46) 



a. Case W >Y 



Using the same decomposition as in the symmetric case (see Figure || and (fn|)), the path- integral (112) becomes 



P?{x,W)=- 



w+r 



sinh ST Jw 

+00 

tlx 



dV m I dx m B r 



(V m ,x m ) (dvoF^ m _ rVm] (W,x\Vo,x m )J \v =v m 



m ( d v F lv m -ry m] ( w > x 'm\Vi x )) W=v m / dF I dlE±(F, I; V m , x' m ) 



-00 p-\-oo 



and thus we obtain in Laplace transform 
p r w+r 



PtfaW) 



w sinh \J p + S 2 T 

sinh y/p + 5 2 (T - (V m - W)) sinh 6(T - (V m - W)) r+ °° 



sinh 2 ST 



sinh yjp + S 2 T 



sinh ST 



dFe~ {F ~ Vm coth sr) 



(A47) 
(A48) 

(A49) 
(A50) 



which leads to the final result (113) given in the text. 



b. Case W < r 



two contributions. The first one corresponding to the case where V m > T is the 
ve to integrate over V m G (T,W + T) (instead of V m G (W, W + T)) and we get 

v+r r ; r 2 



As in the symmetric case, there are two contributions, 
same as ( A4S ) except that now we have to integrate over 

P t {1)± ( P ,W < T) = -^J— [ W+X dV m e* ST ^ P , + 52 e -(Vm-r)(±J+v^cathyOTr) 
* W ' ' sinh 2 ST Jr sinhVP+^r 

* _i_ f~- i F>/-n fir Txrw _• i_ cVti /t; tt7-\\ *»4-qo 



sinh ypT^(r - (K, - w)) sinh^r - (k, - wp) /-+ 00 dFe _ (F _y m)(±5+5cothar ) 

sinh y/p + (5 2 r sinh 5F 7 Vm 



sinh 2 JT | £±g T 4£(t^ + y/pTS 2 coth VpT^r) ) 

V sinh 2 yp+<5 2 r y 

[-2(=F<5 + VpT^ 2 coth ^Tpr) e - w ^ s +V^c°^V^r) 

+2(=F<5 + + S 2 coth x/p+^r) cosh Vp + S 2 W cosh <5W 

-(Vp + ^l + coth 2 y/p + S 2 T) T 25 coth y/p + S 2 T)) sinh a/p+5 2 W cosh 5 W 

+ t(p ^ (T<5 + \/p + S 2 coth ^P + <5 2 r) 2 ) cosh VpT^ 2 sinh (514 7 




(A51) 
(A52) 
(A53) 

(A54) 
(A55) 
(A56) 

(A57) 



-r( T 8 + \/p + 5 2 coth y/p + S 2 T){~2 coth xfp + <5 2 F + ^ ± = ) sinh yJp + 5 2 W sinh <W] (A58) 

<5 sinh y/p + S 2 T 

The second contribution corresponding to the case where V m < T (see Figure (H) can be obtained as in ( p8| ) 
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P^ 2)± (x, w < r)= 



2 



J dlE±(T,z;F,l) 

and thus in Laplace transform 



vb=o 



r2 . /-oo 

* rf^F r y ° F ^] (w ' p|K,) )v n= o/, dz (^' F [tr](^°l W )) y , 

»+oo 



Vb=o Jx 



v=r 



(5 sinh <W sinh(r - W) + 5 2 
sinh 2 <5T sinhr v /p + (5 2 



(A59) 
(A60) 

(A61) 
(A62) 
(A63) 



The final result for W < T is obtained by summing the two contributions (A58) and (A63) and is given in (120) in 
the text. 
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